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Abstract: We analyze the effective 3 dimensional theory previously con- 
structed for the MSSM and multi-Higgs models to determine the regions 
of parameter space in which the electroweak phase transition is sufficiently 
strong for a, B + L asymmetry to survive in the low temperature phase. We 
find that the inclusion of all supersymmetric scalars and all 1-loop corrections 
has the effect of enhancing the strength of the phase transition. Without a 
light stop or extension of the MSSM the phase transition is sufficiently first 
order only if the lightest Higgs mass Mh ^ 70 GeV and tan(3 < 1.75. 
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For electroweak baryogenesis to occur it is necessary that the electroweak 
phase transition be sufficiently strongly ffist order. Otherwise, sphaleron 
transitions after the phase transition wash out any baryon asymmetry which 
may have been produced at the electroweak scale It is now known that 
this requirement is not satisfied in the minimal Standard Model Thus, 
it is of interest to investigate extensions of the minimal Standard Model. 
Many authors have studied the order of the electroweak phase transition 
in the Minimal Supersymmetric Standard Model (MSSM). Most of these 
studies rely on a one- and two-loop finite-temperature effective potential 
analysis of the phase transition |^, ||, |^, || in which stops were expected 
to make the most significant contribution from supersymmetric particles. 
The authors of these studies, in the limit of a large pseudoscalar Higgs mass, 
rriA oo, have identified a region of parameter space for which the transition 
is strong enough. This corresponds to low values of tan(3, and values of the 
soft supersymmetry breaking right stop mass, frifj,^, which are small or even 
negative. In reference |0 the analysis was extended for the full range of 
allowed values of m^. It was found that larger values of rriA are favored. 

A different approach consists of separating the perturbative and non- 
perturbative aspects of the phase transition. This is performed through the 
perturbative construction of effective three dimensional theories, and a sub- 
sequent lattice analysis of the reduced theory ||, ^, ^, For the case in 
which the reduced theory contains a single light Higgs field, characterized 
by a Higgs self-coupling, A3, and an effective 3D gauge coupling, g^, the 
condition for a sufficiently strong first order phase transition becomes 

a;, = ^<0.04, (1) 
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where the quantities A3 and are functions of the various parameters ap- 
pearing in the original 4D theory. 

An analysis of the parameter space for the reduced theory of the Standard 
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Model was performed in 0, 0]. It indicated that for no value of the Higgs 
mass is electroweak baryogenesis possible. This demonstrates that a purely 
perturbative analysis is inadequate, since with that method the electroweak 
phase transition was found to be sufficiently strong for small values of the 
Higgs mass. In |]ll| a 3D theory for the MSSM was constructed, including 
Standard Model particles and additional corrections arising from gauginos, 
higgsinos and all squarks and sleptons. Here we use these results to explore 
the MSSM parameter space in order to determine the regions for which elec- 
troweak baryogenesis may occuifj. This relies on the relation between the 
running parameters in the original 4D theory and physical parameters; this 



is given in |12|. A simplified version of the present analysis was performed 



by others [^3|, |Tj| . However only the contribution from gauge bosons, higgses 
and third generation quarks and squarks to the 3D reduction were included, 
and one-loop corrections to 4D parameters were not fully incorporated. In 
our work all one-loop corrections as well as contributions from all SUSY 
particles have been considered. This allows us to investigate the effect of 
the full complement of supersymmetric particles, in addition to third gener- 
ation squarks, on the strength of the phase transition. We find that these 
effects are important and should not be neglected. We also clarify the dis- 
crepancy between the results of references |jl3| and fl^ . In reference the 
results agreed basically with those found in the perturbative effective poten- 
tial analysis. The most favorable region of parameter space was found to be 
"fTT-h < "^VK (low tanf3), small stop mixing, m^/g < 50 GeV and > 200 
GeV. In addition to this region, reference found another region of pa- 
rameter space in which arbitrary values of tan(3 and a range of values for 
the pseudoscalar Higgs mass, 40 < < 80 GeV, give a sufficiently strong 
phase transition. We comment on the latter region below. 

The ratio Xs/g^ appearing in equation (HD depends on the param- 

^Throughout, we work in the sin^0w = approximation, which was found in ref. ||^ 
to be adequate for the MSM. 



2 



eters in the 4D theory {xc = x{MA,Tno, fi,mi,mg, A,tanf3,Tc)), as well as 
the Standard Model gauge couplings. A is the scalar trilinear soft SUSY 
breaking parameter, taken to be universal; fi is the supersymmetric higgsino 
mass parameter; rrio, mi^irig denote the common squark/slepton mass at the 
SUSY breaking scale 0, the SU{2) gaugino, and gluino mass, respectively. 
Ma is the physical pole mass of the pseudoscalar Higgs, and tan[3 is the ratio 
of the vacuum expectation values of the Higgs fields in the renormalized zero 
temperature theory, x also depends indirectly on the scale Msusy, the scale 
at which the SUSY boundary conditions on the quartic Higgs couplings ap- 
pearing in the Higgs potential are imposed |jl5| and at which mass parameters 
for squarks/sleptons are specified. 

The critical temperature, T^, is defined to be the temperature at which 
there is a direction in field space at the origin of the Higgs potential for which 
the transition to the minimum of the potential in the broken phase can occur 
classically 0. In 3D lattice calculations 0, |I0[ the critical temperature is taken 
to be the temperature of phase coexistence. These two values of temperature 
are generally close but not identical. The actual temperature at which the 
phase transition occurs lies between these two values. We will remark below 
on the circumstances under which there can be a significant difference arising 
from this distinction. 

Throughout our analysis we will concentrate on the regions of parameter 
space which describe an effective theory in which there is a single light scalar 
and thus the bound given by equation (|l|) is valid. However, we mention that 
another possibility is a scenario in which two scalars, e.g. one Higgs and the 
right stop, are both nearly massless at Tc 0. A lattice calculation for this 
extended 3D model is required before that scenario can be investigated with 
^We consider non-universal squark masses where indicated. 

'''We have checked that the difference in the critical temperature from the diagonal- 
ization of the Higgs mass matrix, equation (10) in and from equation (7.9) in 10 is 
extremely small (< .1 GeV) and for our purposes negligible. 
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the present approach. 

As is well known, the MSSM Higgs sector can be parametrized in 
terms of two quantities: tan[3 and the pole mass of the pseudoscalar Higgs 
boson, Ma- These are the most important parameters in determining the 
strength of the phase transition. We consider M4 values between 40-300 GeV 
in order to be compatible with experimental limits and ensure the validity 
of the high-temperature expansion. In the figures we show results for tan(5 
between 1.25 and 13.3. For larger and smaller values the conclusions are 
essentially the same as for the extrema of the range 0. 

In general, the masses of all particles are taken such that the high tem- 
perature expansion is validQ. The experimental constraints we impose on the 
masses are: for stop masses m^^ ^50 GeV, m^^ ^""^t, for first and second 
generation squarks m^. > 200 GeV, sleptons mi > 50 GeV, the gluino mass 
either < 1 GeV or > 150 GeV [0, [l^. In addition, the value of the left soft 
supersymmetry breaking stop mass, m^g, must be such that the contribution 
from stops and sbottoms to the p parameter is not too large [^]. 

Figure |1] shows the value of Xc, for the case of no squark mixing, as a 
function of the pseudoscalar Higgs mass for values of tan(3 ranging from 1.25 
to 13.3. We have fixed the other parameters to be rrio = 50 GeV, mi = 50 
GeV, m-g = ^mi, M^eak = m, Msusy = 10^^ GeV Q. For large values 
of tanj3, there is no value of Ma for which Xc fulfills the condition given by 
equation and Xc varies very little as a function of Ma- However, for low 

^In the dimensional reduction procedure, the explicit dependence on all Yukawa cou- 
plings was kept. For the numerical results presented here only the top Yukawa coupling 
is kept. We have explicitly checked that even for large values of tan/S the bottom Yukawa 
coupling can be neglected. 

^The adequate suppression of non-rcnormalizable terms must also be verified. 

^In our approximation the masses of sleptons and of the first and second generation 
squarks are fixed by nio, rrig, mi, through the renormalization group running; they are 
constant as we vary tanP and Ma- However, due to their dependence on the renormaliza- 
tion group running of the top Yukawa coupling, the left and right stop masses change as 
we move on the curves plotted in figure |l|. 
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values of tan(3 and large enough values of M4, Xc can be small enough for 
the phase transition to be sufficiently first order. The strong dependence of 
Xc on the value of the pseudoscalar Higgs mass, for low tanf3, arises basically 
through the dependence of the quantity A3 in equation on the mixing 



angle 6 (see equation (18) in [|TT1). It is easy to see to lowest order the same 
dependence on Ma arising in finite temperature effective potential analysis 
[^]. This qualitative dependence of the strength of the phase transition on 
tanP and M4 was observed in ||, |l^. Varying the parameters A, /i, nig, 
mi, nio, rrit either increases Xc slightly or has a negligible effect. 

We have also compared the results of our general analysis to those ob- 
tained with the simplifying approximations used in [1T3| , [1^, in which some 
super symmetric particles are neglected. In all cases we have kept the full 
one-loop corrections to the 3D couplings, in contrast to [T^, |14[. Figure || 
shows the variation of Xc with Ma for three different cases and for two val- 
ues of tanf3. The solid line corresponds to our general analysis, as given 
above. The dashed line is the result when only third generation squarks are 
included; the gluino and electroweak gaugino thermal screening contribution 
to the 3D masses of the squarks is excluded. The dotted line corresponds to 
the case in which we include all squarks and sleptons, but ignore all gaugino 
contributions to the three dimensional theory. Although the effect of the 
right stop on the strength of the phase transition is greater than that of any 



other sfermion, the ensemble of sfermions neglected in |]13|, significantly 
strengthens the phase transition. As expected from the work of j^, ^ 
we find that the reduction of the right stop soft supersymmetric breaking 
mass decreases Xc- However, the decrease is less important than it appears 



when only the contribution from third generation squarks |TJ] is included. 
We have also compared the two cases in which only third generation squarks 
were included with and without thermal screening arising from the gluino 
and gaugino. The differences in the values of Xc for this case are negligible. 
For given squark masses at the SUSY scale, the running masses at low 



5 



energies are reduced as the gluino mass decreases. As a result of these lower 
masses the sfermions' favorable impact on the phase transition is increased. 
Thus light gluinos can be helpful in providing a sufficiently strong ew phase 
transition, and the low values of tan(3 required for the phase transition lead 



to chargino masses in an acceptable range in the light gaugino scenario [[Tq] . 

An important point, which has been overlooked in the previous literature, 
is that for some regions of parameter space it may be incorrect to conclude 
from the above analysis that the phase transition is not sufficiently strongly 
first order. As noted previously the actual transition temperature is some- 
what higher than Tc as defined above. For some values of M4 and tanj3, Xc 
depends strongly on the temperature. This happens when the mixing angle 
6, which diagonalizes the 3D Higgs mass matrix at finite temperature, varies 
rapidly for temperatures near the critical temperature. A rapid variation of 
the mixing angle occurs when the temperature is such that the diagonal ele- 
ments of the 3D Higgs mass matrix become nearly degenerate. If the critical 
temperature for the phase transition is close to the value of the temperature 
where this rapid variation occurs, then the value of 6 and consequently of 
Xc are very sensitive to the transition temperature. Note that in general our 
procedure of integrating out the heavier Higgs is not compromised when this 
phenomenon occurs because the eigenvalues remain well-separated. Rather, 
our inability to analyze this region of parameter space arises from our inabil- 
ity (with present techniques) to obtain a sufficiently accurate determination 
of the phase transition temperature. 

The value of the temperature at which this rapid variation occurs depends 
strongly on the value of the pseudoscalar Higgs mass. As M4 increases this 
temperature also increases; for > 100 GeV it is well-separated from the 
phase transition temperature. Moreover, the extent of the variation of Xc is 
less for larger tan(3. Figure ^ shows the dependence of 9 on the temperature 
for Ma = 40 GeV; the solid line corresponds to tan[3 =1.25 and the dashed 
line to tan(3 = 13.3, keeping all other parameters fixed. Figure H is the same 
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for Ma = 300 GeV. Figure || shows Xc as a function of temperature close 
to Tc, for Ma = 40 GeV. A 5 GeV variation in the temperature around Tc 
induces, for Ma = 40 GeV, a change in the value of Xc, Axc ~ .13, while 
for Ma = 300 GeV, Axc ~ 0.005. Thus the possibility of large uncertainty 
in the mixing angle is relevant only for low values of tan(3 and Ma- Since 
Ma ^ 100 GeV is already ruled out experimentally |]l9l for the MSSM in the 
tan(3 region of interest, this possibility cannot enlarge the viable region of 
parameter space in the MSSM and we do not pursue it further. However we 
note that this phenomenon may play a role in the discrepancy between the 
conclusions of ^ H for Ma < 100 GeV. 

We now turn to implications of these constraints for the mass of the 
lightest Higgs. The light curves in figure || indicate contours of constant M/j 
in the tanf5 — Ma plane, using the results given in |]T2| to relate the running 
parameters of the MSSM analysis to the physical parameters and taking 
m^^ ~ 180 GeV, m^-^ ~ 320 GeVQ. Due to approximations in the RG analysis 
and our ignorance with respect to the mass of the stop, one should attach 
a few-GeV uncertainty to these curves. The region of tan(3 — Ma which 
preserves a baryon asymmetry generated at the weak scale is below and to 
the right of the thick line, which corresponds to Xc = 0.04. We conclude that 
unless Mh ^70, electroweak baryogenesis is not viable in the MSSM. 

Given that experimental Higgs limits already nearly exclude such a light 
Higgs, we briefly mention three alternate extensions of the SM which might 
allow electroweak baryogenesis. In order to increase the strength of the phase 
transition one must either increase the (p^T term in Vg/f, or decrease the 
coefficient of the 0^ term, or both. The first possibility is employed in theories 
which are fine-tuned such that one or more scalars in addition to the usual 
Higgs scalar are nearly massless at the phase transition temperature. Terms 
in Veff cubic in arise from bosonic mass-squareds which are proportional 
®This choice of parameters maximizes the region in tan(3 — Ma giving Xc < 0.04. 
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to 0^. Such contributions are enhanced if there is a precise cancelation of 
the thermal contributions to the mass, e.g., from a negative mass-squared 
at the SUSY scale Other scalars such as additional Higgs or sneutrinos 
could in principle be employed to serve a similar purpose, although the stop 
is particularly natural and has 3 color degrees of freedom as well. Non- 
perturbative effects in such scenarios cannot be analyzed without further 
lattice calculations. For the light stop scenario, the 3D theory which must 
be analyzed on the lattice is considerably more complicated than the one 
relevant to the SM and the theories analysed here. The relevant 3D 
lattice calculation for the light stop scenario must include SU(3) as well as 
SU(2) gauge interactions, and one of the scalars couples to both SU(2) and 
SU(3) gauge bosons, as well as to the other scalar. See ref. [jl4|, [Tl[ for the 
full 3D reduction with two light scalars. 

The other strategy to enhance the phase transition is to decrease the 
coefficient of the 0^ term, i.e., A3 in the 3D theory. This may be possible 
either in non-SUSY multi-Higgs doublet models or in the NMSSM (minimal 
SUSY augmented with a gauge singlet Higgs). Non-SUSY multi-Higgs theo- 
ries have some advantage in this regard, since the Higgs potential of the 4D 
theory is not fixed by gauge couplings at the SUSY scale. It is of course also 
constrained by non-observation of a Higgs particle and the requirement that 
the broken-symmetry vacuum is the minimum of the T=0 theory. A dis- 
advantage of non-SUSY theories is the absence of sfermions, which enhance 
the strength of the transition, as we saw above. Finally, the NMSSM has 
sfermions and also more freedom in the Higgs sector, so the lightest Higgs 
in the 4D theory may be acceptably heavy even with a small tan(3, without 



requiring a large Ma 12^, For this theory in particular the 9 dependence 
noted above may prove important in the analysis. 

In conclusion, we have employed existing lattice calculations to analyze 
the electroweak phase transition in the MSSM, including non-perturbative 
as well as perturbative thermal effects. We include all 1-loop corrections 



8 



and integrate out all gauginos and sfermions. Although we find qualitative 



agreement with the results of refs. |13, M, we find that inclusion of all 



sfermions and D-term couplings is important quantitatively, because of their 
large multiplicity. This enhances the strength of the phase transition in the 
relevant regime of parameters. We find that the MSSM provides sufficient 
suppression of sphaleron transitions in the broken phase only for values of 
tan/? < 1.75, unless the right stop soft-SUSY breaking mass is less than 50 
GeV. In the latter case coefficients of non-renormalizable terms become large, 
signaling the onset of the breakdown of our analysis when there are two light 
scalars at the phase transition. In this regime our method does not apply, 
although purely perturbative analysis leads one to expect that the strength 
of the phase transition is enhanced p, 0]. The strength of the phase transi- 
tion increases as the mass of the pseudoscalar Higgs increases. Ma can be 
as low as 100 GeV, for tanfi ~ 1.25, and still give Xc^0.04, even assuming 
universal soft supersymmetry breaking masses at the SUSY breaking scale 
(i.e., without the light stop scenario). The region of parameter space po- 
tentially supporting electroweak baryogenesis requires the lightest physical 
Higgs mass to be < 70 GeV. This is very close to being experimentally 
excluded. We commented on possible alternatives to the MSSM in which 
electroweak baryogenesis could be compatible with experimental constraints 
on the Higgs mass. 

Acknowledgments: We thank M. Shaposhnikov for many useful discus- 
sions. We also thank M. Laine for comments on the manuscript and P. 
Janot for e-mail correspondence. Research supported in part by NSF-PHY- 
94-23002 and Colciencias, Colombia. 

References 



9 



[1] V.A. Kuzmin, V.A. Rubakov, and M.E. Shaposhnikov. Phys. Lett., 
B155:36, 1985. 

[2] K. Kajantie, M. Laine, K. Rummukainen, and M.E. Shaposhnikov. Nucl. 
Phys., 466:189, 1996. 

[3] K. Kajantie, M. Laine, K. Rummukainen, and M.E. Shaposhnikov. Nucl. 
Phys., 458:90, 1996. 

[4] G.F. Giudice. Phys. Rev., D45, 1992. 

[5] J.R. Espinosa, M. Quiros, and F. Zwirner. Phys. Lett, B307:106, 1993. 

[6] A. Brignole, J.R. Espinosa, M. Quiros, and F. Zwirner. Phys. Lett., 
B324:181, 1994. 

[7] M. Carena, M. Quiros, and C.E.M. Wagner. Phys. Lett., B380:81, 1996. 

J.R. Espinosa. Nucl. Phys., 475:273, 1996. 

[9] K. Kajantie, K. Rummukainen, and M.E. Shaposhnikov. Nucl. Phys., 
B407:356, 1993. 

[10] K. Farakos, K. Kajantie, K. Rummukainen, and M.E. Shaposhnikov. 
Nucl. Phys., B425:67, 1994. 



[11] M. Losada. Technical Report ( |hep-ph/ 96052661) , 1996. 



[12] M. Losada. Ph.D. T/ieszs Technical Report ( |hep-ph/9612337D , 1996 



[13] J. Cline and K. Kainulainen. Technical Report ([hep-ph / 9605235|) , 1996. 
[14] M. Laine. Technical Report ( |hep-ph/9605283| ), 1996. 
[15] H. Haber and R. Hempfiing. Phys. Rev., D48:4280, 1993. 
[16] F. Abe et al. Phys. Rev. Lett, 69:3439, 1992. 

10 



[17] G.R. Farrar. Phys. Rev., D51:3904, 1995. 
[18] G.R. Farrar. Phys. Rev. Lett, 76:4114, 1996. 

[19] ALEPH Collaboration. Technical Report EPS-0415, Contribution to the 
International Conference on High energy Physics, Belgium, 1995. 

[20] M. Drees. Int. Journal of Mod. Phys., A14:3633, 1988. 

[21] J. Ellis, J.F. Gunion, H.E. Haber, L. Roszlowski, and F. Zwirner. Phys. 
Rev., D39:844, 1989. 



11 



0.175- 
0.1^ 
0.12$, 



0.075 
0.05 
0.025 



1^0 150 



200 250 300 



Figure 1: Plot of vs. M4 for several different values of tanjS. The solid line 
corresponds to tani3 = 13.3, the dashed line to tanf3 = 1.75, the dashed-dot 
line to tan(3 = 1.5, and the dotted line to tan^ = 1.25. 
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Figure 2: Plot of Xc vs. Ma for tan(3 = 1.5 and 1.75 including all supersym- 
metric particles (solid line), only third generation squarks (dashed lines) and 
all squarks and sleptons but not gaugino corrections (dotted line). 
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Figure 3: Plot of 9 vs. T for Ma = 40 GeV. The solid line corresponds to 
tanp = 1.25, = 63.3 GeV, the dashed line is for tanfi = 13.3, = 75.6 
GeV. 
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Figure 4: Plot of 9 vs.T for Ma = 300 GeV, Tc = 60, 75.6 GeV. The solid 
line corresponds to tanP — 1.25, the dashed hne to tanP — 13.3. 
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Figure 6: Contours of constant Mh = 70,60,45 from top to bottom in the 
tan(5 vs Ma plane. The thick hne gives the constraint for a sufficiently strong 
first order phase transition, = 0.04. 
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